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It is known that in single scalar field inflationary models the standard curvature pertur¬ 
bation C, which is supposedly conserved at superhorizon scales, diverges during reheating at 
times (j) = 0, i.e. when the time derivative of the background inflaton field vanishes. This 
happens because the comoving gauge (p = 0, where p denotes the inflaton perturbation, 
breaks down when (/) = 0. The issue is usually bypassed by averaging out the inflaton os¬ 
cillations but strictly speaking the evolution of C is ill posed mathematically. We solve this 
problem in the free theory by introducing a family of smooth gauges that still eliminates 
the inflaton fluctuation p in the Hamiltonian formalism and gives a well behaved curvature 
perturbation C, which is now rigorously conserved at superhorizon scales. At the linearized 
level, this conserved variable can be used to unambiguously propagate the inflationary per¬ 
turbations from the end of inflation to subsequent epochs. We discuss the implications of 
our results for the inflationary predictions. 


I. INTRODUCTION 

Any successful theory of the early universe must explain the origin of the scale free cosmological 
perturbations seeding the structure in the universe. As a promising candidate, inflation offers a 
natural mechanism that yields a scale free spectrum out of the quantum vacuum fluctuations. The 
well established cosmological perturbation theory studies both the quantum origin and the classical 
evolution of these perturbations (for a comprehensive review see e.g. [1]). An important feature 
that greatly simplifies the evolution of the cosmological perturbations is the existence of a variable, 
the curvature perturbation which is conserved at superhorizon scales and thus directly connects 
the observed fluctuations to the very early times. Indeed, under mild technical assumptions it is 
possible to show that there is always an adiabatic solution to the field equations which is constant 
at large wavelengths [2] (a general proof of the conservation of the curvature perturbation without 
appealing to perturbation theory has been given in [3]). 

In this paper, we consider single scalar held inhationary models where the nearly exponential 
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expansion is followed by the period of reheating during which the inflaton field coherently oscil¬ 
lates about the minimum of its potential. In these models, one repeatedly encounters instants in 
reheating at which p + P = 0, where p and P are the background energy density p = + V 

and the pressure P = — V. As is well known, the conservation of the curvature perturbation 

breaks down exactly when p + P = 0 (see e.g. [1, 2]). Our aim here is to resolve this issue in the 
single scalar field models at the linearized level. 

A very convenient reformulation of the cosmological perturbation theory is given in [4], where 
the curvature perturbation ^ appears as a specific variable in the metric. After fixing the time 
reparametrization invariance by imposing the so called comoving gauge (p = 0 (ip denotes the 
inflaton fluctuation), the scalar degree of freedom is carried out by ("j which is naturally conserved 
at large wavelengths. Indeed, this conservation law can be attributed to a (global) scale invariance 
and it is valid even in the full quantum theory [5, 6] (see also [7]). However, the comoving gauge 
p = 0 breaks down when 0 = 0, i.e. when the time derivative of the background inflaton held 
vanishes, which unfortunately occurs frequently during reheating (note that for a minimally coupled 
scalar held 0 = 0 gives p-|- P = 0). This problem manifests itself as a singularity in the C equation 
of motion, where one sees that a generic smooth initial data diverges in time as 0 ^ 0 and the 
Bunch-Davies mode function is not an exception. Strictly speaking, the evolution of C is ill dehned 
during reheating and it fails to be the demanded conserved quantity. 

In studying the possible effects of reheating on the cosmological perturbations, one unavoidably 
faces the breakdown of C- The standard way of smoothing out the singularities (or as it is sometimes 
referred, the spikes) of 0 is to carry out some form of averaging over inhaton oscillations. For 
example, an important problem is to see whether the oscillations of the inhaton background causes 
a real superhorizon evolution of the cosmological huctuations and 0 is the key variable in answering 
this question. In that case, 0 can be smoothed out by replacing the function 0^ by its time 
average over the inhaton oscillations < 0^ > (see e.g. [8-10]). Note that this averaging avoids 
p + P = 0, since the mean pressure generically vanishes during reheating. In a different setting, 
the same procedure can be used to cure the entropy mode loop corrections to the scalar power 
spectrum during reheating [11]. However, this method is not mathematically well justihed due to 
the substitutions like 1/0^ ^ 1/ <0^>. Of course, it is possible to dehne other variables which are 
well behaved during reheating like the gauge invariant Sasaki-Mukhanov variable or the inhaton 
huctuation held p that carries the scalar mode in the 0 = 0 gauge. The main problem in using 
these variables is that they are not necessarily conserved at large wavelengths and consequently they 
involve nontrivial superhorizon evolutions during reheating. As a result, the details of the reheating 
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or other subsequent stages become important in determining the cosmological observations. 

In this work we take the background inflaton oscillations at reheating seriously, and try to find 
an honestly conserved variable, which is also well defined during reheating. As noted above, and as 
we will discuss below, the main reason for ( becoming an ill defined variable is the breakdown of the 
comoving gauge (p = 0 at times (/> = 0. Thus, the main problem is to find a smooth gauge condition 
that still eliminates the inflaton fluctuation p from the dynamics. The Hamiltonian formalism is 
best suited for this analysis and we indeed find a smooth family of gauge conditions that can be 
used to remove p from the system, leaving a well behaved and conserved curvature perturbation 
C,. Having obtained this smooth variable, we also find transformations that directly relate it to 
the old singular C, and to the Newtonian gravitational potential in the longitudinal gauge, which 
determines the power spectrum of the density fluctuations. As we will discuss, our results put some 
of the standard inflationary predictions on a firmer mathematical basis. 

II. SINGLE SCALAR FIELD MODELS AND COSMOLOGICAL PERTURBATIONS 
IN THE COMOVING GAUGE 

In this section, following [4] we review the derivation of the equations of cosmological perturba¬ 
tions in the single scalar field models. We consider a generic model with the potential U((/>), where 
(j) denotes the inflaton field that is minimally coupled to gravity. In such a model, the background 
field equations can be written as 

(f) 3H(j) + Vff, = 0, 

3H^ = + V, (1) 

H = --(^2 

where = dVjdcj), the background metric is given by 

_j_ ( 2 ) 

H is the Hubble parameter H = dja, dots denote ordinary time derivatives with respect to t and 
we set SttG = 1. For now, the background evolution is completely arbitrary, i.e. we make no 
assumptions like the slow-roll approximation, and any set of fields obeying (1) is allowed for the 
upcoming analysis. Later on we will focus on inflation that is followed by reheating. 

Using the ADM formalism, the action governing the dynamics of this theory can be written as 

S = ^jVh N _ 2y _ ^ (a,,-A b _ ^ , (3) 
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where the metric is parametrized as 

ds^ = -N^dt^ + hijidx^ + N^dt){dx^ + N^dt), (4) 

is the Ricci scalar of hij, Kij = \{hij — DiNj — DjNi), K = Kij and Di is the derivative 
operator of hij. Note that the indices in these expressions are manipulated by the spatial metric 
hij. For a homogeneous inflaton field and the metric (2), one may easily obtain (1) from the 
action (3). 

The cosmological perturbations can be introduced as 

hij — Cl e ^(e'^)ij, 'yij — Tjo la — 0, 

(t) = (t){t) + ip, (5) 

N = l + n, 

N* = 0 + n^ 


where a{t) and (/>(t) are assumed to obey the background equations (1). Note that ^ij is defined 
to be traceless and the determinant of hij is parametrized by Q. At the linearized level, the 
infinitesimal coordinate transformations that are generated by the vector field = {k,k^) yield 
the following gauge transformations 


5C = Hk + \dik\ 

o 

dip = (pk, 

5n = k, 

Sn^ = k' — \:dik, 

5-fij = dikj + djki - ‘^5ij{dmk'^), 


( 6 ) 


where the indices are raised and lowered by the Kronecker delta 5ij. 

Using (5) in (3), one gets the action for cosmological perturbations. In this paper, we only 
consider the free theory and thus it would be enough to determine the quadratic action. The 
equations for the lapse n and the shift n* are algebraic and these can be used to determine n and 
n* in terms of the other fields. On the other hand, the gauge freedom ( 6 ) can be fixed completely 
by imposing 


ip = 0, daij = 0. (7) 

This is called the comoving gauge where the scalar mode is carried out by C and as usual 7 *^ 
represents the gravitational waves. In that case, the lapse and the shift can be solved at the linear 
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order as [4] 


n = 


if’ 


n 


= -d-^di 
2 


9 ^2 _^ 

a?H 


d'^C 


( 8 ) 


Note that the solution for the shift n* involves the Green function for the flat space Laplacian 
9^ and thus it is non-local in the position space. Nevertheless, ( 8 ) can easily be understood in 
momentum space. As discussed in [4], to obtain the action up to cubic order in fluctuations it is 
enough to determine n and n* at the linear order. One may now use ( 8 ) back in the action (3) 
to eliminate n and n* completely, leaving a Lagrangian for the three physical degrees of freedom 
carried out by C, and 7 *^. A straightforward calculation shows that the quadratic actions decouple 
5 ( 2 ) = 5 ® _l_ 5 "^^) ^ where 



It is important to note that no slow-roll approximation is used in deriving these actions, i.e. they 
are valid provided that the background field equations ( 1 ) hold. 

The full set of equations must be invariant under the transformation a —)■ Aa that corresponds to 
a constant scaling of the spatial coordinates x. From (5), where the perturbations are introduced, 
the scaling a —)■ Aa can be interpreted as a shift of C as C C + This shows that the 

constant C configuration, which is generated by this shift from C = 0 ; should solve the C equation 
of motion to all orders in the perturbation theory [4]. This special property leads to the appreciated 
superhorizon conservation law. 


III. THE BREAKDOWN OF C DURING REHEATING 


One sees from the action (9) that the evolution of C becomes singular at (j) = 0. Indeed, the 
field equation of C can be found as 


c + 


3H + 


2 ( 1 ) 


2H 

IT 


C- —dX = 0, 


( 11 ) 


where the second term in the square brackets involving 1/0 is problematic: Given an initial generic 
data set C(^ 0 ) x) and ({to, x), the solution developed by ( 11 ) is singular where C, necessarily diverges^ 


^ One also sees from the action (9) that the propagator of ^ blows up in quantum theory at </ = 0. 



6 



FIG. 1. During reheating the inflaton oscillates about the minimum of its potential and at the turning 
points A and B, one has (p = 0 causing a problem for the C evolution. In time, the amplitude of oscillations 
decreases i.e. A and B approaches to the origin. 


in time at (p = 0. For the class of inflationary models we are considering, one frequently encounters 
instants for which (p = 0, see Fig. 1. 

The singular behavior emerging from (11) can nicely be illustrated for the superhorizon modes. 
Introducing the standard Fourier space components (pk, the two solutions of (11) for negligible 
wavenumber k can be found as 


Ck-C 


( 0 ) 

k 


+ Cfc 


dt' 


a^{t')(p‘^{t') ’ 


( 12 ) 


where and are two constants, which are fixed by the vacuum chosen and the time of horizon 
crossing Since cp is oscillating, one has (p ^ t) near (/> = 0, and thus the integral in the second 
solution diverges^ like f dx/x^. One may suggest to discard the second solution since 

it is decaying during inflation and becomes completely negligible for the modes of cosmological 
interest (the Wronskian condition, which is required for the canonical commutation relations to 
hold, implies — Ck^^*<^k = h and thus it is not possible to set Ck = 0). However, this decaying 

solution naturally appears in the loop effects of entropy modes during reheating [11]. Moreover, 
the modes leaving the horizon just before inflation and reentering the horizon in reheating might 
have important physical effects (for instance, they may cause primordial black holes to form, see 
e.g. [12]) and for those modes the decaying solution is not completely negligible. 

On the other hand, the “constant” mode is not safe either. The first order derivative correction 


^ While the integral f dx/x can be cured by taking the principal value, this method does not work for J dxjx^ . 
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to the superhorizon mode can be obtained from (11) as [13] 


Among other things, this correction determines the semiclassicality of the superhorizon modes [14], 
and we see that it once more blows up as the first integral passes through the singularity of the 
integrand at = 0 . 

In general, it is easy to prove that any smooth generic initial data set, which is evolved by (11), 
would yield a solution that diverges like C oc l/4> near 0 = 0. To see this, one may introduce a new 
variable by (^ = — (0/i7)0. An easy calculation then shows that the field equation of (p becomes 
completely regular even at 0 = 0. Therefore, ( must blow up like 1/0 as 0 ^ 0. Indeed, (p is 
nothing but the inflaton fluctuation, which is an entirely well defined variable. Using the gauge 
freedom ( 6 ), an arbitrary configuration {‘p,C) can be transformed to ( 0 ,C) or (<^, 0 ) by choosing 
k = or k = —(/H, respectively. As a result, the gauge transformation from (0,0) to (v?, 0) 

gives the relation ip = —(0/77)0. These comments also explain the main reason for the breakdown 
of 0. Namely, to impose the comoving gauge (^ = 0 one must take k = —pjcj) in (6), which gives a 
diverging gauge parameter at 0 = 0. Therefore, the comoving gauge is not an allowed condition if 
0 = 0 . 

IV. A FAMILY OF SMOOTH GAUGES 

In this section, we discuss how one can eliminate the inflaton perturbation p from the dynamics 
by imposing smooth gauge fixing conditions. Whether such smooth conditions exist and whether 
they can be used to remove p from the system are not evident in the Lagrangian formulation, 
where the comoving gauge (34) looks like the only possible option. Therefore, we switch to the 
Hamiltonian formalism where the gauge invariance is connected to the existence of constraints in 
the phase space (see e.g. [15] for a similar Hamiltonian analysis). This would give a better geometric 
picture about the dynamical evolution of the system and it may suggest better alternatives to the 
comoving gauge. 

Since we would like to work out gauge fixing in the Hamiltonian formalism, we keep all 11 
variables 0, p, n and n* in the Lagrangian. Note that is traceless and thus it contains 
only 5 independent degrees of freedom. Plugging (5) into (3) without imposing any conditions, the 
quadratic action of all fluctuations (recall that in this paper we only consider the free theory) can 



Ck 


AO) 


1-r 


dt'- 


H^{t') 


02 ( 0 ') 


a3(0)02(p) 


dt''a{f 


H^{t" 



be obtained as 


/ II 1 -1 

- -^a{dajk)idajk) + -^a{daik){dj'yjk) - Sa^C^ + adiCdiC + 

--adi(pdi(p - -a V(p + -aCdidjjij - 3a <j)Cp (14) 

-^a^diUj'jij + 2a^Cdini + ^a^{dinj){dinj) - ^a^{dinif - a^^Uidip 

-o?n —^didjjij - — 2V sP —+ 12HC — 4:Hdini —a^Vn^, 

2 [ a^ J 

where = dVjdci), all index manipulations are carried out by the Kronecker delta 5ij (like n* = n*) 

and an over-line indicates that the function is evaluated in the background solution. The canonical 

momenta conjugate to 'jij, (, p, n and n% which are respectively denoted by 11*^, P^, Pn and 

Pi, can be found from this quadratic action as 

In Inf 2 

7b “ ~ -^^ijdknk , 

P(^ = —6a^ C “ 5a^4>p + 2a^dini + 6Ha^ n, 

P^p = p — a^^n, (15) 

Pn = ^, 

P, = 0. 

One can then apply the standard Legendre transformation to obtain the Hamiltonian^ of the system 
as H = J d^xP: 

P P g 0, di'yjkdi'yjk ~ ci dj'yjkdi'jik CLdiQdi(^ I2a^ 2a^ 

\ _ 3-1 1 

+ lja?V^^p^ - + -adipdip- -aCdidjjij -n4>-n*4>i, (16) 

where 

= 2djllji + —diP^, (17) 

<1> = - a didjjij — 2ad'^Q + (pp — ^Pp, — HPq. (18) 

Since not all components of ^ij are independent due to trace-freeness, it is convenient to introduce 
a basis in the space of 3 x 3, traceless, symmetric matrices, which we denote by T7, a, 6 = 1,2,5. 
One may normalize these matrices to satisfy 


rjia rjib _ ^ab 

rpa rpa _ ^ 

ij ^ mn 2 


^im^jn “ 1 “ ^in^jm 


^ While H is the Hamiltonian, Ti denotes the respective density. With some abuse of terminology, we also refer Ti 
as the Hamiltonian. 
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We expand jij and Ily as 7 ^^ = jaT^j and liij = liaT^, and treat 7 a and Ha as independent 
variables. 

The canonically conjugate variables obey the standard Poisson bracket relations; 


= 5^{x - y), 


{^p{t,x),P^{t,y)] = - y), 

1 , 

{lij{t,x),Yimn{t,y)} =-5 (x-y) 




jn 


+ 6 , 


in^jm 




ij 'Jmn 


( 20 ) 


where in getting the last relation we have used the properties of the matrices given in (19). The 
constancy of the primary constraints in (15) under the Hamiltonian flow, i.e. Pn = 0 and Pi = 0, 
implies 


$. = 0, ^> = 0. (21) 

As usual, n and n* become non-dynamical Lagrange multipliers. On the other hand, the first class 
constraints and ‘h, which are are related to the gauge invariance ( 6 ), obey 

= 0 . ( 22 ) 


A direct calculation shows that the constraints are preserved under the Hamiltonian flow; 


dt 

dt 






'm 


= 0 . 


(23) 


This is a highly nontrivial consistency check of the expressions given above. One must note that 
the constraint has an explicit time dependence through the background helds a{t) and 
which are assumed to obey the equations of motion (1). For a given vector field = {k,P) one 
may define & = f d?x{k^ + which generates the gauge transformations in ( 6 ) 

5C = {Q,C} = Hk + U,k\ 

(jc/? = {0,99} = #, (24) 

Slij = { 0 , 77 } = dikj + djki - ‘^ 5 ij{dmk'^). 


The transformations of the conjugate momenta can be found as 


SP^ = { 0 ,PJ = -2ad‘^k, 


6P^ = 

SUij = 


{0,P<^} = a^^k 


{ 0 ,n,,} 



didjk 


-diid'^k 
3 ^ 


( 25 ) 
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Since we treat n and n* as non-dynamical Lagrange multipliers, their gauge transformations can 
be imposed from the invariance of the action S = f (pq — H) under (24) and (25) (see e.g. [16]), 
which gives the corresponding equations in (6), i.e. 

5n = k, 

SrC' = F — \dik. (26) 

One may check that the gauge transformations (24), (25) and (26) are consistent with (15). 

After obtaining the Hamiltonian and the first class constraints, we now proceed with the gauge 
fixing. As is well known, a function that has a non-zero Poisson bracket with a first class constraint 
defines an allowed gauge for the corresponding invariance. After finding such a suitable gauge 
condition, one may prefer to work either in the full phase space by replacing the standard Poisson 
brackets with the Dirac brackets [17], which are defined by the solutions of the Lagrange multipliers, 
or rather “solve the constraints” to obtain a reduced phase space and a reduced Hamiltonian [18]. 
The equivalence of these two procedures, when there is no explicit time dependence, has been 
proved in [18]. In the appendix A, we generalize that result in a simplified but time dependent 
setting, which is suitable for our discussion. 

The gauge freedom corresponding to can conveniently be fixed by imposing 

Gi = = 0, (27) 


which obeys {^i,Gj} 0. This condition completely fixes the spatial diffeomorphisms generated 

by H. To preserve Gi under the Hamiltonian flow, i.e. to have dGi/dt = {Gj, 77} = 0, the Lagrange 
multiplier n* must be set to 

To make the constrained phase space defined by the conditions $* = Gj = 0 more obvious, one 
may use the decomposition 

m = iW + + didj-i - 

Hi, = + 25(,nJ) + didjH - Uijd^Ii, (29) 

where 'yj'/' and HJ]-^ are transverse-traceless tensors, and 'yf and are transverse vector fields. 

Ij Ij 'll 

The equations $* = G* = 0 imply 



( 30 ) 
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The motion is now confined in this (partially) constrained phase space defined by (30). To obtain 
the reduced Hamiltonian Ti! generating the dynamics in this subspace, one may use (30) directly 
in (16) (see the appendix A). Not surprisingly, the gauge condition (27) decouples the dynamics 
of the tensor and the scalar sectors from each other giving l-i' = + T-is-, where and Hs are 

the respective Hamiltonians. We find that the Hamiltonian of the tensor modes becomes 

^ a dajk, (3i) 

which corresponds to the standard action for the gravitational waves (10). 

On the other hand, the scalar sector involving the fields C, and (p has the Hamiltonian 

Tis = -adiCdiC - + ^adipdip + ]^o?V (32) 

where the first class constraint is given by 

^ = <j)ip — ^ ^ (33) 

As a consistency check, one may verify that d^/dt = {^,Hs} + d^/dt = 0. Note that for 
convenience we rescale the constraint in (32) by as compared to (16). Time reparametrizations 
produced by the vector fields of the form = (/c,0) are generated by 0 = J d?xa?k^. Gauge 
fixing can be done in various ways and as we discuss in the appendix B it is possible to rederive 
the previously known results from (32) by imposing (/? = 0 or ^ = 0 gauges. 

We would like to find a smooth gauge condition that completely eliminates the inflaton fluctua¬ 
tion from the system, leaving C, as the main scalar mode. Since (33) determines a linear combination 
of p and Pip, one should impose a complementary condition that can be used to solve for p and 
Pp. We set 

G = c{t)p + d{t)Pp, = 0, (34) 

where c{t) and d{t) are two arbitrary functions. To have a unique solution for p and Pp from (33) 
and (34), one needs 

det ^ / 0 (35) 

c d 

and without loss of any generality we impose 

c + = 1, 

O'* 


(36) 
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which normalizes the above determinant. During reheating when 0 = 0 one has 0 7 ^ 0, and vice 
versa. Consequently, it is possible to choose completely regular functions c{t) and d(t) obeying 
(36) (see (65) for an explicit example). Since {^{t, x),G{t,y)} = 5^{x — y), the function G defines 
a viable one parameter family of smooth gauge conditions.^ The constancy of G, i.e. dG/dt = 
{G, Hs} + dG/dt = 0 determines the lapse uniquely as 

n = [<«%♦ - + ^ + - di) - [^P( + 0 ay] - (37) 

As we will see below, this expression greatly simplifies when expressed in the configuration space. 
It is straightforward to eliminate if and by using (33) and (34): 



As we discuss in the appendix A, the reduced Hamiltonian for the curvature perturbation^ 0 can 
be obtained by both plugging these solutions into the Hamiltonian (32) and adding the extra term 
emerging from a time dependent canonical transformation, which is derived in (AlO). By comparing 
(33) and (34) with the corresponding equations (Al) and (A3), one sees that the functions c(t) and 
d{t) are identical and by comparing (33) with (Al) one determines the function K that appears in 
(AlO) as 

K = ^- • ( 39 ) 

Using (38), (39) and (AlO), the reduced Hamiltonian for 0 can be found as 

%(■ = -adiCdiC + ^d^diKdiK + ^dKd‘^C + AT^, (40) 

Z li Z 

where the time dependent function y, is given by 

(b \ 3 1 r 1 

11 = + ^ ^ dc-dc . (41) 

H Z a'^ L J 

Note that (40) contains spatial derivatives of the momentum variable P(^ through (39), which is 
unconventional. 

^ Note that given d{t), the other function c(t) can be determined from (36). Below, we also reinterpret (34) in terms 
of the configuration space variables. 

® At the linearized level, Q still determines the Ricci scalar of the induced metric hij on the constant-t hypersurfaces 
since di'yij = 0. 
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Let Cc be the curvature perturbation in the comoving gauge that has the action (9). One can 
see that the following time-dependent linear canonical map 

H‘^d 


Cc = C- 


(f)a? 
2a 


Pc + 


Pc. =PC + 

transforms the Hamiltonian of Cc to (40), which proves the equivalence of the dynamics. 
The quadratic action corresponding to the Hamiltonian (40) can be obtained as 


(42) 


c(2) _ 


1 


where 


-M 


M = 




a(j4 


diCdiC, 




H- -^d 


-1 


(43) 


(44) 


In the Fourier space, M must be evaluated by the substitution 8“^ —)■ —k'^, which fixes the notation 
in (43).® As we give an explicit example below, it is possible to choose c and d so that ^ > 0 and 
M becomes strictly positive. Although the action is nonlocal in the position space, it involves only 
two time derivatives and the propagation for is well defined once the initial conditions C{to,x) 
and ({to,x) are given at some time to. Eq. (43) generalizes the action (9) of the comoving gauge. 
Indeed, setting d = 0 in (36) gives c = a^/(j) and (43) reduces to (9). 

In obtaining the action (43) from (40), a timelike integration by parts is applied. Due to the 
surface term thrown away in (43), the original canonical momentum P(^ is given by 


Pc = 


- —d^C = M 
SC H 


C - ^d^c 


-|a^ 


(45) 


This relation will be important in the next section. 
The equations of motion can be found as 

Ft 



d^ 

= -yM 


L V « yJ 






(46) 


As long as c and d are chosen smoothly, p. and M become well behaved and one encounters no 
singularity in the C equation during reheating. Not surprisingly, the constant configuration 


C = Co (47) 

solves (46). As pointed out above, this is because of the fact that C is related to the global scaling of 
the metric function a —)■ \a. Therefore, with the gauge condition (34) the curvature perturbation 
C becomes a well defined conserved variable during reheating. 

® Namely, the Green function must be understood to be sandwiched between the two terms in (43). 



14 


V. IMPLICATIONS 

From the expression for given in (15), one may see that the gauge condition (34) can be 
written as c(p + a? d[Lp — (l)n] = 0. Using this equation, it is easy to see that given an arbitrary set 
of perturbations ((^, n) one can apply the gauge transformation (6) with the parameter 


k = -X 


ap + a? d [(f — 


(48) 


so that (34) is satisfied (recall that the functions c and d are normalized to obey (36)). As a result, 
the gauge conditions introduced in the previous section can be written as 


dnij = 0, C(f + d 


ip — (pn 


= 0 . 


(49) 


We refer (49) as the cd-gauge, which completely fixes the infinitesimal diffeomorphism invariance 
in the free theory. In this gauge, the basic variables are {'jijX) and the other perturbations p, n 
and n* can be expressed in terms of C- Indeed, using the momentum in (45), the lapse n can 
be fixed from (37) as 


C 


Similarly, (28) determines n* as 






i-*4s\ 


(pd 


(50) 


Pa2 


diC 


(51) 


and (38) implies 


p = —M 

nO 


c-4b\ 


(52) 


Note that only the derivatives of C appear in these expressions for n, rd and p. 

As emphasized above, the constant field C = Co solves the equations of motion. Neglecting the 
spatial derivatives^ in (46), the second “superhorizon” solution can also be obtained so that 


Ck = C1+C2 / dt /c ~ 0. 


(53) 


As long as c oc and d is chosen independent of a, one can see from (41) that p becomes 
independent of a. Thus, the two solutions in (53) can be referred as the constant and the decaying 
modes, and unlike the solution given in (12) these are completely regular during reheating. 


It is possible to see that (46) admits a derivative expansion for ^ 7 ^ 0. 
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It is possible to transform the curvature perturbation ( in the cd-gauge to the curvature pertur¬ 
bation in the comoving gauge, which has been denoted by Cc- Since the comoving gauge is defined 
by the condition (p = 0, one can apply a coordinate change to the perturbations in the cd-gauge 
that sets (p = 0. The corresponding gauge parameter k in (6) can be found from (52), and applying 
the coordinate change gives® 


Cc 


H^d 


M 



(54) 


where C is the curvature perturbation in the cd-gauge. Similarly, starting with the comoving gauge 
with (p = 0 and the lapse given in (8), one can apply a coordinate transformation to generate new 
fields p and n obeying (36). Finding the corresponding parameter and changing the curvature 
perturbation yield 


C — Cc + d^Cc- 


(55) 


Note that (54) and (55) are inverse of each other provided that C and Cc obey their respective 
equations of motion. We check that the lapse n and the shift n*, which are determined by the 
respective curvature perturbations, also transform accordingly under these transformations. 

Since C is smooth, (54) shows that Cc diverges like 1/C> as i;A ^ 0, as it should be. On the other 
hand, defining a new variable v = <pC,c, (55) becomes C = dv + cv/a^, which shows that C is smooth 
since v is well behaved. 

The power spectrum of Cc can be obtained by quantizing the action (9) on the inflationary back¬ 
ground. In the slow-roll approximation with the Bunch-Davies vacuum, the constant superhorizon 
mode, which is denoted by in (12), becomes (see e.g. [1]) 


Co. -cf = 


1 


(56) 


^{tk) 

where tk is the time for horizon crossing: k/a{tk) = H{tk) (the scalar power spectrum is defined 
by ICfc^^P)- Assuming further that the solutions in (53) give the constant and the decaying modes 
respectively,® (55) shows that to a very good approximation the constant superhorizon modes are 
the same in the comoving and the cd-gauges: 


Cfc = Ccfe = cf, 


/c ~ 0. 


(57) 


® Eq. (54) can also be obtained from (42) and (45). 

® As noted above, this depends on the behavior of c and d. For the functions given in (65) below, this identification 


is true. 
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Namely, during inflation at superhorizon scales the initial C spectrum is identical to Cc spectrum. 
As can be seen from (55), changing the function d only affects the decaying solution for the super¬ 
horizon modes. Consequently, our new variable C, has the standard inflationary power spectrum 
and it safely propagates this initial scale-free spectrum beyond reheating. 

An important variable that has a direct physical meaning is the curvature perturbation in the 
longitudinal gauge, or the Newtonian gravitational potential^® O, which determines the density 
perturbations.^^ In the longitudinal gauge, the metric takes the following form 

ds^ = —(1 — 2(^)dt^ + (1 + ‘2<t))a^d3p‘ (58) 


i.e. in terms of the perturbations introduced in (5) one has 

C = O, n = —(h, n* = 0. 


(59) 


It is possible to transform C, or C,c to <!> by finding the suitable parameter k that sets n* = 0 in the 
transformation (6), which gives 


O = 


f M 

2 a 


2H O' 


d-‘^C - gC 


cd — gauge 
comoving gauge. 


(60) 


Eq. (60) shows that the superhorizon mode is fixed by the leading order derivative correction 
of the constant superhorizon curvature perturbation and an easy calculation gives 




^ / <■(* 


^( 0 ) 


/c ~ 0. 


(61) 


H{t'Y 

Note that the Newtonian potential O is well defined at all epochs. 

Let us hnally consider an explicit example for the cd-gauge that would illustrate some of the 
points above. We first specify the background evolution for a reasonably large set of generic models. 
During inflation, the background dynamics is controlled by the standard slow-roll parameters 
defined by 

1 yl Vs. 


e = 


V = 


(62) 


2 1 ^ 2 ’ ■' Y 

We take a massive inflaton with mass m, therefore the potential in reheating can be taken as 
V ~ ^rn^cj)^. Assuming further that H, which is generically true in many models, the scalar 
held equation in (1) can be solved approximately as 


(j){t) ~ <I>sin(mt), (63) 

To be precise, the Newtonian potential arising from the linearization of general relativity around flat space is — O. 
The relation between O and 6p can easily be obtained from the perturbed Einstein’s equations. 
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where the slowly changing amplitude obeys <l> + 3i7<h/2 ~ 0. In that case the background equations 
yield the following approximate solution 

$ ~ j^o, (64) 

which is valid during reheating. 

A convenient choice for the functions c and d is 



m? (j) 


c = 


d = 


<P 


(65) 


+ (j)^ TvP'cjP' + (fP' 

which satisfies the normalization condition (36). During the exponential expansion, these functions 
are determined by the slow-roll parameters as 


c ~ - 


d ~ 


rj — e 




On the other hand, (64) can be used to obtain the form of the functions in reheating 


c ~ 


cos(mt) 




sin(mt) 


( 66 ) 


(67) 


’ m^<I> ’ 

which are completely well behaved. From (65), the second super horizon solution given in (53) can 
be seen to be decaying like 1/aP during inflation. 

Another suitable choice is (recall that we are using the Planck units SttG = 1) 




+ (^2 ■ 


d = ^ 


'2 + 02 


( 68 ) 


Using (63) and (68) in (41), we numerically check for a various phenomenologically interesting set 
of parameters that one has a strictly positive function fi > 0 during reheating. 

In this model, one can use (61) to relate the Newtonian gravitational potential to the initial 
inflationary power spectrum. Using (63) in (61) and remembering that we have assumed m S> 77, 
the highly oscillatory integral can be approximated by treating the slowly changing functions as 
constants that gives 


^ (69) 

which is the standard relation between and in a matter dominated universe. This calculation 
proves that there cannot be any amplification of the superhorizon density perturbations 6pk if the 
infiaton is massive since is constant. This is a highly nontrivial result that has been shown in [8] 
by carefully analyzing the equations of motion for perturbations. We see that the same conclusion 
can straightforwardly be reached by using the conserved variable (). 
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VI. CONCLUSIONS 


The conservation of the snperhorizon curvature perturbation is crucial in directly relating the 
late time cosmological observables like the CMB temperature perturbations to the quantum me¬ 
chanical vacuum fluctuations in the very early universe during inflation. The general belief is that 
due to the conservation law the details of the cosmic evolution after inflation does not matter for 
the late time properties of the cosmological perturbations. However, it is not much appreciated in 
the literature that the standard curvature perturbation in the comoving gauge becomes singular 
during reheating if the inflaton oscillates about the minimum of its potential. 

In this paper, we try to solve this problem by introducing smooth alternative gauges, which are 
both regular in reheating and still eliminate the inflaton perturbation from the dynamics. As we 
have shown, this can be achieved in the Hamiltonian formalism where the coordinates and momenta 
are treated as independent variables. Although it obeys an unconventional equation of motion, the 
new curvature perturbation becomes a smooth variable in reheating whose superhorizon mode is 
rigorously conserved. 

It turns out that the new gauge can be chosen in such a way that the constant snperhorizon mode 
in inflation is equal to the corresponding mode in the comoving gauge. In other words, only the sub¬ 
leading decaying pieces differ between gauges. Therefore, the new curvature perturbation has the 
standard inflationary power spectrum and using its conservation law it is possible to propagate the 
initial scale free spectrum beyond reheating. For instance, an important relation is (69) that relates 
the Newtonian potential and hence the density perturbation to the initial inflationary spectrum. In 
early works (see e.g. [9]), (69) has been obtained from the conservation of the (singular) curvature 
perturbation by approximating the coherent inflaton oscillations by a matter dominated universe 
on the average. Here, we obtain this result in a rigorous way without referring to the singular 
variable. 

In this work, we have focused on the reheating era and assumed that the perturbations during 
inflation can be treated in the standard way. However, there has been some concerns about the 
behavior of the cosmological perturbations in the limit (/> —)■ 0 as the nearly exponential expansion 
approaches to the exact de Sitter space, since the standard power spectrum (56) diverges (see 
e.g. [19]). For the moment, our findings cannot be directly applied to analyze this issue since we 
assume that (j) ^ 0 when <p = 0 and vice versa. Nevertheless, it would be interesting to work out 
the necessary modifications in an attempt to clarify (/>—)■ 0 limit of inflation. 
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Appendix A: Gauge fixing and phase space reduction in the constrained Hamiltonian systems 

In this appendix, we consider a classical mechanical system whose dynamics is determined by 
a Hamiltonian and a first class constraint. We assume that the constraint is at least linear in a 
coordinate and the conjugate momentum, and otherwise arbitrary. We study the evolution after 
gauge fixing in two equivalent ways; either by using the Dirac’s method of the constrained systems 
[17] or by reducing the dynamics on the constrained phase space. Our analysis is a simplified 
version of [18], however we generalize that discussion by taking both the Hamiltonian and the 
constraint to have explicit time dependencies, which is suitable for cosmological applications. We 
show that the two methods are equivalent, i.e. they yield identical evolutions. 

Let us denote the Hamiltonian as H{q,p,qi,pi]t). The first class constraint is assumed to have 
the form 


^ = a{t)q + b{t)p-K{qi,pi;t), (Al) 

where K is an arbitrary function and a{t) and b{t) denote time dependent external parameters. 
Being a first class constraint, <I> must obey d^/dt = {^,iL} + d^/dt = 0. The total Hamiltonian 
that generates the most general motion can be written as 

Ht = H + A$, (A2) 

where A is a Lagrange multiplier that is completely arbitrary at the moment. 

Due to the existence of the constraint, the evolution of the system is actually confined in a 
codimension one sub-manifold defined by <I> = 0. Besides, <I> generates a gauge transformation and 
one may fix that freedom to reduce the motion in a codimension two manifold, which can be viewed 
as a reduced phase space. The gauge freedom can be fixed by introducing a function G that obeys 
{G, 7 ^ 0. In that case the codimension two surface in which the motion is confined is given by 

$ = G = 0. For our gauge condition we take 

G = c{t)q + d{t)p = 0 (A3) 

and without loss of any generality impose a(t)d{t) — b{t)c{t) = 1 so that {<!>, G} = 1. As mentioned 
above, to study the dynamics of the system, one may either follow Dirac, which amounts to 
determine the Lagrange multiplier A from the condition dG/dt = 0, or solve for the conditions 
(Al) and (A3) to find the reduced codimension two phase space and the corresponding reduced 
Hamiltonian. 



20 


In the first case, dG/dt = {G, Ht} + dG/dt = 0 fixes A as 
1 


A = — 


{G,H} + ^ 


{G,^} 

and the equations of motion become 


dH ^ dK 
9 * — - 1 “ —) 

upi ^Pi 


GiTJ 

=-dt)^ + dit)-^ - c{t)q - dit)p, 


dH ^ dK 
Pi = “ o- 

uQi UQi 


(A4) 


(AS) 


The evolutions of q and p are determined by the conditions ‘h = G = 0, which read q = 
d{t)K{qi,pi]t) and p = -c{t)K{qi,pi]t). 

Let us now try to obtain a gauge fixed Hamiltonian H{qi,Pi;t) that describes the dynam¬ 
ics of the “true physical degrees of freedom” g* and pi. Of course, the flow generated by 
H{qi,pi;t) must be identical to (AS). Since = G = 0 implies q = d{t)K{qi,pi;t) and 
p = —c{t)K{qi,pi',t), one may first try a direct substitution in the Hamiltonian (A2), which 
gives H{d{t)K{qi,pi]t), —c{t)K{qi,pi]t),qi,pi]t). But it is easy to see that the equations generated 
by this Hamiltonian misses the terms involving c{t) and d{t) in (AS) through (A4). Here, the 
explicit time dependence of the gauge condition causes a problem. 

To resolve this issue, one may first apply the following canonical transformation 


Q = a{t)q + b{t)p, 

P = c{t)q + d{t)p, (A6) 


which can be generated by the function 

= + (A7) 

where one has p = dFi/dq and P = —dFi/dQ. Under this transformation, the Hamiltonian must 
also be shifted as 

dPt 

H'{Q,P,qi,pi;t) = H{q{Q, P),p{Q, P),qi,pi) + —, (A8) 

where the substitutions in Fi must be done after the partial time derivative is evaluated. In these 
new canonical variables {Q, P,qi,pi), the constraint and the gauge fixing condition G become 

Q = 

P = 0. (A9) 


In these equations the new canonical coordinates, which are going to be eliminated, are not multi¬ 
plied by time dependent factors and one may obtain the reduced Hamiltonian on the constrained 
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surface by simply using (A9) in (A8), which yields 

= H {d{t)K{qi,pi;t), -c{t)K{qi,pi;t),qi,pi;t) + ^ d{t)c{t) - c{t)d{t) K{qi,pi;tf. 

(AlO) 

It is easy to check that the equations generated by this Hamiltonian are identical to (A5). Eq. 
(AlO) is the basic result of this appendix, which is used in the main text above. 

Appendix B: Rederiving the equations in the previously known gauges 


In this section, we rederive the well known equations corresponding to the standard gauges 
p = 0 and C = 0 from (32). Let us start by imposing p = 0 which corresponds to the comoving 
gauge discussed in section II. To eliminate the inflaton field completely, one may use the constraint 


to solve for as 




(Bl) 


As we discuss in the appendix A, the reduced Hamiltonian can be obtained by setting p = 0 and 
using the solution for in (32) that yield 

1 r- _ .,2 a12 


= -ad.CdiC + [HP^ + 2ad^C] . 


(B2) 


It is straightforward to show that the corresponding Lagrangian is equal to (9). Moreover, one can 
check that the solutions for the Lapse n, which can be obtained from {(^, Hs} = 0, and the shift 
n* given in (28) are identical to (8). 

Similarly, it is also possible to impose C = 0 and use (33) to solve for as 


H 


Setting C = 0 and using this solution in (32) give 

1 n a ^ ^ 

+ Y 


t7 _ ^ 

<P4> 2^ H 




After a Legendre transformation, one then obtains 


^( 2 ) = 


n 


- ^dipdip - Y W + 


1 


, 3 -2 ^ 2 


(B3) 


(B4) 


(B5) 


a- ■ H 

which is the standard action for the inflaton fluctuation, see e.g. [4]. It can be further checked that 
the shift n* and the lapse n, which are obtained from (28) and from the condition {^, Hs} = 0, 
also match the usual expressions given in [4]. Consequently, we show that it is possible to obtain 
the well known equations for the cosmological perturbations in (^ = 0 and C = 0 gauges from (32). 
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